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Abstract 


A new method for the calculation of shell model intrinsic density matrices, de¬ 
fined as two-particle density matrices integrated over centrum of mass position 
vector of two last particles and accompanied with isospin variables, has been 
developed. Produced intrinsic density matrices are completely antisymmetric, 
translational invariant and do not employ a group theoretical classification of 
antisymmetric states. They are devoted for exact realistic density matrix ex¬ 
pansion within the framework of Reduced Hamiltonian method. The procedures 
on a base of precise arithmetic for intrinsic density matrices calculation that in¬ 
volve any numerical diagonalization and orthogonalization were developed and 
implemented in computer code. 


I. INTRODUCTION 

Last time discovering of a exotic nuclei emphasize the dehciency of usual nuclei description 
methods based on central held approximation. The reason relay upon week exotic nuclei binding 
energy due to oversaturation of neutrons or protons. To describe such the systems large-space 
shell model expansions have to be considered, since the ordinary shell model is unable to manage 
deep nucleon-nucleon correlations. However shell model functions are implicitly redundant 
with many times encountering a numerous of the same intrinsic states. The reason is that 
excited state shell model wave function is representable as linear combination of the products 
of intrinsic wave function and centrum of mass (c.m.) wave function and only the hrst term in 
the series corresponding to the c.m. in ground state is needful. Every remaining term contains 
the one of excited intrinsic wave functions that was encountered in the hrst term of expansion 
of less excited state shell model wave function. The hrst attempt to solve this problem was the 
translationally invariant shell model preliminary invoked to antisymmetryze the wave- 

function depending on intrinsic coordinates, leads only to sophisticated exercising in group 
theory which do not end in available results. 

We propose a simple and ehective procedures comprising simplicity of ordinary shell model 
and requirement of translational invariancy for wave-function calculation. The method relies 
on the number of statements simplifying the problem under consideration. 

At hrst, despite the harmonic oscillator potential do not depend on spin and isospin co¬ 
ordinates, we refuse of orbital, spin and isospin spaces separation in model wave-functions. 
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Hence we loose an advantage of using the group theoretical methods for many-particle anti- 
symmetrical states classihcation. Nevertheless, on the other hand, rejection of precise labelling 
of model wave-functions by quantum numbers of higher order groups, such as seniority, leads 
to simplification of the very antisymmetrization procedure. In this approach many-particle an- 
tisymmetrical states are characterized only by well dehned set of quantum numbers: number of 
oscillator quanta E, angular momentum J, parity U, isospin T and only one additional integer 
quantum number E = 1,2,.. .r necessary for unambiguous enumeration of the states. Here 
r is rank of the corresponding antisymmetrization operator matrix. Developed on this basis 
computational procedures allow unrestricted conhguration space to be taken into account @,§]. 

Second simplihcation refer to the ’spurious’ states. Description of atomic nuclei can be car¬ 
ried out within the framework of ordinary shell model using intrinsic wave-function expansion in 
terms of shell model functions. The coefficients of this expansion can be obtained diagonalizing 
c.m. Hamiltonian in the shell model basis. Intrinsic motion of nucleons would be represented 
by subspace of c.m. Hamiltonian eigenvectors corresponding to minimal eigenvalue equal |. 

At last, simplihcation of the A-particle system description is using two-particle density 
matrices instead of wave-functions. Since all operators of observables including intrinsic Hamil¬ 
tonian and root-mean-square (r.m.s.) radius operators are symmetric ones, their expectation 
values calculations do not require the complete wave-function to be involved. Moreover two- 
particle density matrix is still redundant in its variables. We can use more simple quantity as 
density matrix integrated over c.m. position vector of two last particles, as so-called intracule 
0. In case of light nuclei intracule have to describe A-particle system with not excited c.m. 
and must be accompanied with isospin variables @]. Such intracule we call intrinsic density 
matrix. Goal is that using intrinsic density matrices we have to deal only with 7-1-7 orbital, 
spin and isospin variables, instead of wave-functions containing 5A — 3 variables. 


II. DEVELOPMENT OF THE METHOD 


In the shell model it is assumed that free nucleons are moving in self-consistent central 
held. It can be approximated by an isotropic harmonic oscillator potential. The many-particle 
Hamiltonian used in shell model is 
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Here m denotes the nucleon mass, r* the i-th nucleon radius vector, a; the angular frequency 
and h{ri) the single-particle hamiltonian. Usually in the shell model only one eigenfunction of 
many-particle Hamiltonian (H) is taken into account. Hence the consideration is restricted by 
only one conhguration (so-called ground conhguration) characterized by minimal total oscillator 
quantum number. The consistent consideration could be achieved only on complete basis 
of eigenfunctions of Hamiltonian (|2.1|) . To construct any antisymmetrical eigenfunction of 
many-particle H (|2.1|) characterized with a good quantum numbers it is appropriate to dehne 
eigenfunctions of single-particle h{fi) in a j-j coupled representation as 


'^eljtrrijmti,^) .ReZ (^){Tj(7) ® tti/2 (cr) 


( 2 . 2 ) 


Here Rei{r) is a radial function, a spherical harmonics, Q:i/ 2 ms(<?) ^ spin-1/2 function 

in spin space and tti/ 2 mt('r) a spin-1/2 function in isospin space. The single-particle variables 
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are Xj = riffiTi (a set of the corresponding radius-vector, spin and isospin variables). The e, / 
and j are the principal, orbital and total angular momentum quantum numbers, mi, and 
mj are the magnetic projection quantum numbers of orbital, spin and total angular momentum 
respectively. The is projection of isospin dehned such that mt = -1-1/2 corresponds to a 
neutron state and mt = — 1/2 to a proton state. A vector coupling of the angle and spin 
functions to form a state of good total angular momentum is denoted by {... ® 

Intrinsic properties of a nucleus can be described only by intrinsic wave function. It depends 
on [A — 1) intrinsic (Jacobi) coordinates and spin-isospin variables. In general eigenfunctions 
of Hamiltonian m could be represented as linear combination of products of intrinsic wave 
functions and c.m. functions when all c.m. excitations are taken into account. Let us determine 
the system of orthonormalized Jacobi variables according to the 2A —1 vertices Jacobi tree 
presented in Fig. 1. The upper vertices correspond to the single-particle variables xi.. .xa- 
The remaining A — 1 vertices (situated below the hrst ones) correspond to the Jacobi variables 
which could be of three types: 


• hrstly if vertex associated with Jacobi variable is not connected directly with any single¬ 
particle variable, the Jacobi variable is taken to be equal the corresponding Jacobi coor¬ 
dinate (as associated with = ,^i) 

• secondly that connected with only one of single-particle vertices is associated with Jacobi 
variable which is composed of Jacobi coordinate and directly reached single-particle spin- 
isospin variables (as ^3 = ^ 3 UA- 2 fH- 2 ) 

• hnally to that connected with two of single-particle vertices will be prescribed Jacobi 
variable which is composed of Jacobi coordinate and two sets of spin-isospin variables 
coming from directly connected single-particle vertices (as ^2 = ^ 2 ^a-iTa-i^aTa)- 


Orthogonal transformation to Jacobi coordinates is 

1 




/ PaQa 
Pa + Qa 




E o-f E 


r,- 


Qa 




, a = 1... A — 1 


(2.3) 




Here Pa is the number of single-particle vertices which could be reached while moving from 
the a-th vertex upwards along the left edge; {pa} is its numbers manifold and {ga}-the 
same for the right edge. The ^0 is proportional to the nucleus c.m. coordinate. For the vertex 
connected with two last single-particle vertices we obtain for instance ^2 = — fA). 
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Fig.l Jacobi tree definition 


In the Jacobi coordinates the Hamiltonian (|2.1|) takes the form 
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( 2 . 4 ) 


It implies that the Hamiltonian could be divided into two parts: c.m. Hamiltonian 
and Hamiltonian iJjntr representing the intrinsic motion of the nucleons in the system. The 
set of h(,^o)’s eigenfunctions TeZm('Co) consist of products of corresponding radial functions and 
spherical harmonics. Whereas h(,^a)’s single-particle eigenfunctions could be of three different 
types depending on number of single-particle variables directly connected with a-th vertex (if 
any): 


• hrstly there are functions depending on Jacobi coordinates which have no direct connec¬ 
tion with any single-particle coordinate (as ^i). Such the functions are products of radial 
function and spherical harmonic 

• secondly there are functions of the form (p.2|) for Jacobi coordinates directly connected 
with one single-particle coordinate (e.g. ,^ 3 ). 

• thirdly there will be functions depending on a Jacobi coordinate and two sets of spin- 
isospin variables (e.g. ^ 2 )- They can be represented as product of the vector coupled 
orbital and the vector coupled two last particles in the list spin functions and vector 
coupled isospin functions of corresponding particles. 
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'>Pelj-KtmjmM 2 ) = 8) {ai/2(<?A-l) ® 01/2 (^a) {ai/2(TA-l )®aV 2(?^)},„ (2.5) 

Here orbital function '0eZm(G) is product of the radial function and spherical harmonic. 
The intrinsic wavefunctions ^ErjnTMjMriii ■ ■ ■ ^a-i) are eigenfunctions of Tfintr.- Since the 
antisymmetrization procedure is rather cumbersome in Jacobi coordinates it is convenient to 
introduce the expansion of the product of c.m. ground state function and intrinsic wave-function 
in terms of shell model functions 


^ErJUTMjMxi.^^ * * * l) ^ ^ ^EKAJIJTMjMt^^^ * * * ^a) A;OO^r (2.6) 

KA 

The coefficients of this expansion r can be obtained diagonalizing c.m. Hamiltonian 

in the shell model basis. The summation in this formula runs over all configurations K 
and additional quantum number A (in spirit of the T). Here zeros indicate the c.m. ground 
state: principal and orbital angular momentum quantum numbers. The intrinsic wave-function 
depends on the 5H — 3 intrinsic Jacobi variables where stands for the together 

with one of the mentioned three types of the sets of spin and isospin variables. 

Due to the identity of the nucleons it is convenient to restrict consideration only by two last 
particles in the list. General form of an expression of two-particle intrinsic density matrix is 
chosen to be 

^'^EEjnTMjMri^l^ ^2 ■ ■ ■ ^A-l)'^*E'r'jnTMjMT(^^^ ^2 ■ ■ ■ 'Ca-i) 

where integration sign denotes integration over the continuous variables of the relative 
motion of H — 2 particles and summation over discrete variables of the corresponding nu¬ 
cleons according to the Jacobi tree in use. Here and in the following the notations as 
[J,T] = (2J -|- 1)(2T -|- 1) is shortcut of number of states with corresponding angular mo¬ 
mentum and isospin. The most convenient way to construct an antisymmetric wave function 
is the method of fractional parentage coefficients (CFP) developed in @]. 


'^EEjnTMjMri^i • • • ^A-i) = E [(ErjUT)-elj7it\{ErjUT] 

(EEJET) (2.8) 

eljirt ^ ' 

^ EE. ® JETMjMt 

Here double bar indicates the grandparent state and eljTrt characterize the separated sub¬ 
system. The summation in this formula spans subspaces of grandparent {EF JUT) and two 
particle e/jvrf states which satisfy the necessary selection rules and the energy, momentum and 
parity conservation conditions. The expansion coefficients have to be intrinsic CFP. By means 
of (p.8|) it is possible to get the expression of the intrinsic density matrix 


QErjnFET'jnT^^^^^,^^^ E E [{ErjnT)-el3Tit\\ErjnT] 

^i^'i'j^^jETjnr) 

x[(Erjnne'i'iEt\\E'r'jnT] y V4««.,.,„.(?2 )Ctj.i,„,„.(6) 

VJ') mj,mt 


(2.9) 


5 


















The intrinsic density matrix in harmonic oscillator representation is 


E [(ErjnTy,dj„t\\ErjnT][{ErjnT)Wt3^t\\E'r'jnT] 

{ErjWf) 

This matrix could be dehned avoiding the appearance of intrinsic CFP. The key is rela¬ 
tion (|2.6|) between intrinsic and shell model wave functions. The usual shell model fractional 
parentage expansion could be used to separate out two last particles in the list 


EKAjnTMjMTi^l ■ ■ ■ ^a) — _ 

^ {{EKAjnT); {{elj)A-i, {elj)A)rr'\\EKAjnT) 


X 


(EKAjnT) 

{elj)A-i,{elj)A,J"T'' 

{^ {EKAjnrj ^^^ • • -^^- 2 ) ® ^iielj)A-i,{elj)A)J"T'' 


{xa-1,Xa) 


JETMjMt 


( 2 . 11 ) 


Here inverted commas refer to the separated two-particle subsystem. For the sake of con¬ 
venience of the notations we shall use comma to denote the nucleons variables with respect to 
which a wave-function is antisymmetric. The coefficients of introduced expansion (generalized 
CFP (GCFP) as dehned by Levinson |^) enable to express the antisymmetric shell model wave 
function of A nucleons in the form of linear combinations of products of antisymmetric func¬ 
tions of H — 2 nucleons and antisymmetric two-particle wave functions. The separation of the 
two nucleons from the initial conhguration can be accomplished in all possible ways consistent 
with required triangular relations, thus giving rise to transformation matrix describing trans¬ 
formation between different momentum coupling schemes. In case when two nucleons are taken 
from different shells the GCFP can be expressed in terms of single shell one-particle CFP and 
corresponding transformation matrix 


{(EKAjnT)-, ((eij)4_i, (elj)AJ"r'\\EKAjnT) = “ 

x(((JiTi ... JrTr ... JpTp . . . JkTk)J T, {j rtr,jptp)J”T”) JT\ 

|(JlTi . . . {JrTr,jrtr)JrTr ■ ■ ■ {JpTp, jptp) JpTp . . . JT) 


Here single bar over the quantum numbers indicates the parent state, subscript r refers to 
the r-th shell in the conhguration and superscript is number of particles contained in the 

k 

r-th shell. The integer number rii, where sum runs over all shells standing to the 

i=r-\-l 

right from the r-th shell. The {{elj)^'^~^{AJT)^-, {elj)r\\{elj)^^{AJT)r) denotes the one-particle 
CFP of r-th shell. When two nucleons are separated from the same shell the GCFP dehnition 
contains two-particle CFP: 


{{EKAJUT); ((e/j)A-i, {elj)A)rr'\\EKAjnT) = 

f riririr — 1) 


{{el3)r-\AJT)^- {el3)lJ"T"\\{el3)r{AJT\) (2.13) 

x{((JiTi... « ... JfcTfc) JT, J''T'')JT\{JyTy ... (X^, rr')JrTr ... JT) 


\ A{A - 1) 
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The simplest way of CFP matrix calculation is proposed in |^. The method is based 
on the observation that the spectral decomposition of antisymmetrization operator matrix 
is not uniquely dehned. The best choice is to set the upper triangle of CFP matrix equal 
to zero. This method was implemented in the program code and all s-d shells CFP were 
calculated. Corresponding two-particle CFP were obtained following the well-known Redmond 
formula. Transformation matrix describing momentum recoupling was calculated using its 
direct representation by sum of Clebsh-Gordan coefficients. 

The intrinsic Jacobi variable ^2 could be introduced by expressing antisymmetric two- 
particle wave-functions ^({eij)A-i,(eij)A)J"T"{.XA-i,XA) in terms of single-particle eigenfunc¬ 
tions (p.2|) and coupling them in proper order for well-known Talmi-Moshinsky transformation 
to apply. At hrst the antisymmetric wave functions ^([eij)A-i,(eij)A)J"T"{.XA-i,XA) should be 
expressed as a linear combination of not antisymmetrized coupled momentum wave-functions. 
Here and below we will indicate such functions with semicolon as separation mark for the 
variables in the list. In case when the nucleons are taken from the same shell we have 


^{(elj)A-i,{elj)A)J''T"{xA-l,XA) “ 2 


J"+T" 


^{{elj)A-A{elj)A)J''T"{xA-l] Xa) ( 2 - 14 ) 


Here we suppress for brevity magnetic quantum numbers. When the two nucleons are taken 
from different shells linear combination is of the form 


^{(elj)A-i,{f^lj)A)J''T"{.XA-l,XA) = [^{(elj)A-i;{elj)A)J"T"{xA-l]XA) 

-(-iyA-i+3A+l-J -T 


(2.15) 


We reveal the complete representation of the coupled momentum wave functions by the 
coupled orbital-spin and isospin functions 


^{(elj)A-i\{elj)A)J''T''M'jM!^{xA-l'i Xa) — {ai/ 2 ('rA-l) ®ai/ 2 ('y 4 )} 

X |{0e^_p^_i(rA-l) ® ai/2{aA-l)}.^_^ {(l>eAlAirA) ® ai/2(aA)} 


I T"M^ 


(2.16) 


J J»M" 


Here orbital function is dehned as (t)eimif) = Reiij)Yim{f)■ 

Second point concerns transition to the appropriate Jacobi coordinates. This transformation 
from single-particle variables xa-i and xa to Jacobi variable ^2 could be accomplished by 
assistance of Jacobi coordinates with nonpositive indices ^j. The Inters are chosen to be 
proportional to the c.m. of corresponding subsystems. 
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1 

yJPa + Qa 


j&{Pay\qa} 


1...A-1 


(2.17) 


For example .^2 is accompanied with ^ 
c.m. of the two last particles. When a - 
index coincide with ,^0 dehned in (|2.3| ). 
and down the Jacobi tree is dehned 


_i = ■^{ta-i + '^a) which is simply proportional to the 
= 1 corresponding Jacobi coordinate with nonpositive 
This orthogonal transformation allowing to climb up 







e-M- 

+ 





(2.18) 
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Here da = Valdai is Jacobi coordinate associated with vertex which could be reached 
while moving from the a-th vertex upwards along the left edge; is the same for the right edge. 
In case when any of the upper vertices is bound to single-particle variables the corresponding 
single-particle coordinate is taken. 

To apply this transformation it is necessary to interchange the order of coupling. It will be 
the so-called L-S coupling scheme. 

® ® {(leAlAi'^A) ® ^1/2 (< 7 ^) | = 

= ^{{Q-A-li^n)3A-li {lAi^/‘i)jA)J''\{{lA-lilA)L, (l/ 2 ,l/ 2 )s)J") (2.19) 

Ls 

X {{(t)eA-VA-i{^A-l) ® (t>eAlAi.^A)} ^ ® {«l/2(<?A-l) ® ai/2(an) } 

Here L is the orbital momentum of the relative motion of the two last nucleons and s is the 
corresponding spin. Now can be used the Talmi - Moshinsky - Smirnov transformation 


LM 


SA-VA-li^A-l) 'i ^ 

-- ((cOa-i, (cOa : J^|(eO-i,e/: L)i {'0(ez)_j(f_i) 

{el)-i,el 


( 2 . 20 ) 


LM 


here sum is restricted under energy {ca-i + = e_i -t- e) and parity conservation. The 

Talmi - Moshinsky - Smirnov coefficients {{el)A-h [d) a '■ T|(e/)_i,e/ : following |1l0|,pd| 

are the elements of the matrix for the transition between the oscillator functions depending on 
coordinates which are related with orthogonal transformation (|2.18|) . Finally we must return 
to the momentum coupling defined in the i>eij 77 tmjmt{^ 2 ) functions 


J"M" 


(eO-lK-l) ® V’d(f2)}j^ ® {ai/2('?/l-l) 0ai/2('?/i)}J 

3 

X V’(eZ)_i(f-i) ® {^/’ez(6) ® \ai/2{dA-i) ® ai/2(^n)} J . 


( 2 . 21 ) 


^^ 3 ) J"M'J 


Now we can express functions from the left side of the expression (p.l9 ) in terms of the 
functions from the right side of the expression (|2.21|) 


{0eA-iZA-i(^A-l) ® ai/2{dA-l)] ® {(leAui'^A) ® ai/2(^A)} 


' ) Jl'M'j 


= H {{{el])A-ii{el])A)J''\{{el)-iielsj)J”) 

{el)-i,elsj 

X W(eZ)_i(^-l) <8) (V'ezfe) <8) \ai/2{dA-l) ® ai/2(^A)}J 


( 2 . 22 ) 




The coefficients of introduced expansion represents the collection of the transformation ma¬ 
trices from (|2.21|) and ( p.l9|) , and the Talmi - Moshinsky coefficients from the expansion (|2.2(]|) 


L 

x{{{Ia-Ii 1/2)JA-1, {Iai l/2)iA)J''\{{lA-lilA)L, (l/2,l/2)s)J") 

x((e/)A-i, {el)A : T|(e/)_i,e/ : L)i 


(2.23) 
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Corresponding transformation matrices can be represented in terms of standard vector co¬ 
efficients: 6-j and 9-j, thus we are led to the expression 


' I A- 

X Y^Lmd)A-l, {el) a) : L|(e/)_i,e/ : L)i 


1 


/_i I L 
s J" 


J 


1 2 
2 

L s r 


(2.24) 


It should be noted that presented coefficients weighted only not antisymmetrical functions 
as is stressed by semicolon. Taking into account the isospin part of the two-particle function 
we get the final form of the coefficients for the transition from the antisymmetrical two-particle 
shell model functions to the sought 'ipeijiTtmjm,t{^ 2 ) vector coupled with the introduced functions 


{{{elj)A-i, {dj)A)rT"\{{el).,,eljnt)rT") = 


\T" 


-1) 




(l + ^0 


hep)jv-i,Cb)iv )) 

x{{{elj)A-i,{elj)A)r\{{d)-i,elsj)r) 


(2.25) 


This transformation allows us to re-express the two-particle shell model function expansion 
in the form 


4> 


((elj)A-i,{elj)A)J"T"M'jM^{xA-l, Xa) — XI {“l/ 2 ( 7(4-1) ® 01/2(7(4)} 

{el)-i,elsj 


T"M” 


X |V’(eZ)-i(f-l) ® {^eK 6 ) ® (01/2(071-1) ( 8 ) 01/2(0^) 

X {{{elj)A-i, (e/j)A) J"T"|((e/)-i, elj7Tt)rr') 


(2.26) 


*7 J J"M'J 


Finally we obtain the eigenfunctions of A-particles Hamiltonian (|2.1|) with separated two 
groups of nucleons: hrst, containing two last particles which are described by antisymmetri¬ 
cal, coupled function, depending on Jacobi coordinates, and second, containing remaining set 
of particles, which are described by corresponding eigenfunction of H — 2-particles Hamilto¬ 
nian depending on single-particle coordinates. That functions we call the A-particle oscillator 
functions with singled out dependence on intrinsic coordinates of two last particles in the list. 

EKAjnTMjMriXl ■ ■ ■ Xa) = 


= E E E {{EKAjnT);{{el).,,eljTTt)J"T"\\EKAjnT) 

(EKAJIIT) FO-i eljitt,J"T" (2.27) 


X 


(EKAJET) 


(Xi . . . XA-2) < 8 ) {V’(eZ)_i (f-l) < 8 ) i^elj nt (&)} 


J"T" 


JETMjMt 


The momentum coupling scheme between the functions can be arbitrary since density matrix 
is invariant with respect to orthogonal their transformation. Taking advantage of this will save 
us for one more recoupling looked for coupling scheme with 'il^eijnt{^ 2 ) coupled at end. 

The coefficients of the (p.27|) functions expansion by introduced ones are composed of GCFP 
and coefficients (^.25|) 
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{{EKAJHT)-{{el)_i,el]nt)J"T"\\EKAjnT) = 

= Y. ^KAjW)-{{elj)A-u{elj)A)J"T"\\EKAjnT) ( 2 . 28 ) 

{elj)A-i,{elj)A 

x(((e/j)A-i,(e/j)A)J"T"|((eO-i,e/j7rt)J"T") 

This set of coefficients may satisfy the orthonormalization condition 

^ ^ ^ {{EKAjnTy,{{el)_,,eljnt)J"T"\\EKAjnT) 

(EKAjnT) hO-l eljnt,J"T" (2.29) 

x{{EKAjnf); {{el)-i,elj7rt)rT”\\EK'A'jnT) = 6 ka,k'A' 


Now we are able to get the hnal expression of intrinsic density matrix. Preliminarily it 
should be noted that due to orthogonality of transformations ( |2.8| ) and ( |2.18 ) it may hold 
equivalency for integration over two sets of variables 


d^od^id^s ... d^A-i = dxi... dxA-id^-i 


(2.30) 


Here the integration is taken in sense of the definition (|2.7|) . Inserting in the density ma¬ 
trix (El the linear combination of shell model functions with fixed c.m. state (0) and 
integrating after substitution of variables according (p.30 ) we get the sought intrinsic density 
matrix expression 


Q 


ErjnT,E'r'jnT 


^EJET E'JET 

-^Qji 


(^ 2 ,^ 2 )— X! X! ®^Aoo,r"'K'z\';oo,r' 

ele'l'jnt KA,K'A' 


X ^ {{EKAjnTy{{el)-i,eljTrt)J"T"\\EKAjnT) 


(EKAJET) 


(2.31) 


X ({ERA JRT ); ( (e/) _ 1 , el'jut) J"T" \ \ E'K'A' JHT) 

[J'^ mj,mt 


M) 


Comparing the expressions at the products of the functions V’eq 7 rtmjmt(^ 2 )t/’eq’ 7 rtm,mt(^ 2 ) 
the intrinsic density matrix expressions (|2.9| ) and (|2.31|) we are led to the final result 


E {{EKAjnTy,ndu,dj 7 ,t)j"T"\\ErjnT) 


(EKAJET) 


(2.32) 


X ({ERA JHT ); ( (e/) _ 1 , e'l'jnt) J"T" \ \ E'E' JRT) 


Presented coefficients are connected with outlined above coefficients (|2.27|) and (|2.6|) ac¬ 
cording to the relation 


{{ERAjnT)-{{el).i,el]Tit)J"T"\\ErjnT) = 

= YUERAJET); ((e/)_i, djnt)J"T"\\ERAjnT) a^^A%,r 

KA 

That coefficients project out the space of A-particle oscillator functions with singled out 
dependence on intrinsic coordinates of two last particles into intrinsic functions subspace with 
not excited c.m.. 
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The intrinsic density matrix in harmonic oscillator representation satisfy a usual normal¬ 
ization relation 


elj nt 


EFjnT.ET'jnT 

eljirt^eljirt 


= 5 


ERET' 


(2.34) 


The presented expressions enable to obtain the full set of intrinsic density matrices describing 
A-nucleon system in isospin formalism after calculating the coefficients ( ^.331) weighting linear 
combination of that eigenfunctions of A-particles Hamiltonian (^) for which ones c.m. is in 
hxed state. 


III. COMPUTATIONAL RESULTS 

Since momentum recoupling, antisymmetrization and transformation to Jacobi coordinates 
are orthogonal transformations a precise arithmetics could be applied. Instead of calculations 
with real numbers, which are connected with serious numerical instabilities, calculations were 
performed with numbers represented in the form n/{my/k), where n, m and k are integers. 

On the basis of precise arithmetics were developed computational procedures of a number 
of well-known transformations coefficients and implemented in a computer code. As follows: 
one-particle single shell CFP [^, two-particle single shell CFP according to Redmond expres¬ 
sion, Clebsh-Gordan coefficients, 6-j and 9-j vector coupling coefficients, momentum recoupling 
matrices, GCFP, Talmi coefficients. 

Another set of computational procedures was developed for enumeration of the antisym¬ 
metric states by means of combinatorial calculations. 

Outlined above computational procedures were used for general formalism, presented in 
the section Number 2, to implement in a computer code. For illustration let us propose the 
calculations of 6, 7 and 8 nucleon systems with minimal oscillator energy compatible with the 
Pauli exclusion principle were performed. The number of calculated intrinsic GFP and intrinsic 
density matrices for 6, 7 and 8 nucleon systems are: 255 and 41, 1345 and 66, 5021 and 138 
respectively. 

As an example we consider the case A = 6. In this case minimal total oscillator quantum 
number is ~ 2. Hence 6 nucleon system can be in the three ground state conhgurations 






/ 1 \ 4 / 1 \ 2 


Ki = 

(ooi) (lli) : 

JT = 01, 10, 

K2 = 

(00l)*(lll) (ll|) : 

JT = 10, 11, 20, 21, 

K^ = 


JT = 01, 10, 21, 30. 


Here (e/j)” denotes the single shell. Displayed sequence of them taken in ascending quantum 
number order stands for conhguration. According to the number of configurations with the 
same JT values, situated at the right, we have to obtain: one three-dimensional matrix with 
JT = 10, two two-dimensional matrices with JT = 01,21 and three one-dimensional matrices 
with JT = 11, 20, 30. All shells and conhgurations are unambiguously characterized by their 
JT values, thus additional quantum numbers are not necessary. The same specihcation scheme 
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by JT values will be valid and for coefficients of expansion of intrinsic wave-function by shell 
model functions. According to the Elliott and Skyrme theorem shell model states characterized 
by contains the nuclear c.m. in its ground state, hence o^'^oor = 1- ^his case the 

coefficients of expansion of intrinsic wave-fnnction by shell model fnnctions are the CFP of 
intrinsic fnnction. 

Let us consider the intrinsic CFP characterized by JT = 30. This is the simplest case since 
they can be generated only in danghter conhguration (^^l) • Computational results 

are presented in Table 1. 

Table 1 

Intrinsic CFP of 6 nncleon system: 

Emin = 2, K= (ooi)" (llf)' , A = 1, JT = 30 


KA 

JT 

(e/)_i 

elsjt 

state 

J"T" 

CFP 

(OO0'O1(11|)'"3O 

31 

00 

00001 

'^0 

01 

1 

vT 

10 

20 


00110 


10 

_ ]_ 


30 





_ 1 
v45 


40 





3 






vT05 

TooiT^WniT^ 

V 2^ 22 12^ 22 

10 

11 

00110 


20 

1 

2j5 


11 

00 

11121 

"P2 

"Pl 

21 

3 

2vT5 


20 


11010 

10 

_ ]_ 

J30 



11 

00110 



1 




2vT5 






20 

1 

J30 


21 

00 

mil 

"Pl 

11 

1 

2j5 



11 

00001 

"So 


1 

do 



00 

11121 

"P2 

21 

1 

do 

(ooD'oo 

00 

00 

22130 


30 

1 

J30 



22 

00110 



1 

a/1) 


Composition of grandparent conhguration KAis clearly seen from displayed notations where 
every single shell is specihed by its total JT valnes. Proposed in this paper set of qnantum 
nnmbers TJIITMt for enumeration of antisymmetrical states of A-particle system can be in 
natural way prolongated and for only two particles. The corresponding set will be Ijntmt. 
Here orbital momentum I plays the role of addition quantum number as T, since in the case 
of two particles / is well dehned quantum number. So we follow the more usual spectroscopic 
notation to denote the state of the subsystem containing two separated particles. To 

avoid overloading the Table 1 we show only the valnes of the indices which in turn were not 
repeated in the preceding row. 

The obtained intrinsic CFP were used for intrinsic density matrices calculations according 
the dehnition (2.31). As an example we present intrinsic density matrix characterized by 
E = = 2, JT = 10 when the snbsystem of two separated particles is in state S'! 
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-1 


(3.1) 


W23»(®Si) 


44 1 -1 

135 135 27^10 

44 1 

135 27v^ 

17 
54 


Here we simplify notations of (2.31) and display in superscripts only E and JT. ^Si) 

is symmetric matrix, thus only upper triangle is shown. It rows and the same columns are 
labelled by intrinsic CFP left hand bracet indices. Namely rows labels are: 
hrst - ((00|)^00, (ll|)^10)l0, 

second - (((oOi)^OO, (ll|) (llf) |i)l0, 

third - ((oOi)^OO, (ll|)^10)l0. 

Here notations of Table 1 are used except that parentheses represents the momentum cou¬ 
pling order. Using the ordinary shell model notations for intrinsic CFP is justihed by their 
clarity, origin of representation space construction procedures and preservation of space dimen¬ 
sion under transformation to Jacobi coordinates. Last means that spaces labelled by E and 
KA have equal dimensions. 


Table 2 

Diagonal elements of 6 nucleon system density matrices 


JT 

10 

30 

01 

20 

21 

11 

Conhg. 

term 

(O' 

(0 (i) 

iir 

(1)' 

(O' 

(O' 

(0 (0 

(0(0 

(O' 

(0 (0 


44 

44 

17 

1 

3 

3 

1 

3 

3 

3 

135 

135 

54 

3 

1Q 

1Q 

3 

10 

10 

10 

" 5 ; 

1 

8in 

8 

405 

1 

81 







3 

3 

3 

3 

14 

29 

3 

29 

14 

3 

10 

10 

10 

10 

45 

90 

10 

90 

45 

10 

i q ' 



1 

90 

1 

45 




"Pi 

13 

13 

19 

1 

1 

1 

1 

1 

1 

1 

270 

270 

270 

30 

30 

30 

30 

30 

30 

30 

"Po 

1 

1 



1 


1 

1 


13 

in 

20 



10 


20 

20 


180 

"Pi 

1 

1 

1 

1 

11 

13 

1 

47 

11 

17 

5 

8 

20 

20 

45 

180 

8 

360 

180 

120 

"P2 


1 

1 

1 


1 

1 

17 

17 

11 


8 

4 

4 


4 

8 

120 

60 

72 






1 

45 

1 

9Q 



2 

81 

1 

lfi2 

1 

405 






"D2 




1 

30 





1 

30 






The empty places in the Table 2 means that there are any intrinsic density matrices of such 
kind. Here for the brevity in the conhguration notations we don’t display the hrst closed shell 

and suppressed e and 1 in the higher single shells notations. 

Presented results imply that even in approximation we can make same predictions 

concerning effects invoked by tensor character of nuclear forces. The display of variety effects 
arise from so-called mixing of states and have the same origin as in deuteron case. According to 
the reduced Hamiltonian method ^ presence of bound states of A-particle system is explained 
by existence of two-particle bound states of the so-called reduced Hamiltonian. Direct numerical 
integration of reduced Hamiltonian equations with some realistic nuclear interaction potential 
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could supply corresponding eigenvalues and percent of eigenstates mixing. For A = 6 that are 
only ^Sq and ^Si — ^Di. The admixture of ^Di term may consist about 30% . This admixture 
must be exhausted by corresponding terms with higher orbital momentum in expansion of 
realistic density matrix by intrinsic ones. The less part of admixture is exhausted in lover 
energy approximations the more remain to the higher ones. Hence less pronounced contribution 
of Di term in approximation forces increasing influence of admixed terms arising from 

higher energy approximations. Increasing contribution of higher energy approximations in turn 
must increase the calculated nuclear interaction radius. 

As an example we consider the nuclei ^Li and which ground states are characterized 
by JT = 10 and JT = 01 correspondingly. As can be seen from Table 2 the admixture of 
^Di term is absent for ^He whereas it is not the case for ^Li. Hence influence of higher energy 
approximation must be more pronounced in the ^He case and as consequence it must have the 
larger interaction radius, whitch is in full correspondence with experimental results. 


IV. CONCLUSIONS 

The only quantities really needed for calculation of identical particle systems are two- 
particle density matrixes. The developed formalism is particulary pointed to description of 
energy spectrum formation mechanism and irregularities of root-mean-square radii (Halo effect) 
in light exotic nuclei. The method consistently lines the principles of antisymmetrization and 
translational invariancy and is implicitly based on Reduced Hamiltonian method. The proposed 
density matrices are very suitable for exact realistic density matrix expansion since it enable 
do not involve the realistic spectator functions, devoted for description of the remaining A —2 
particle subsystem, into consideration. It is well-known that such oscillator many particle 
functions form slowly convergent series and thus the as large as possible set of states have to be 
used. This way the procedures of group theory for choosing the ’best state’ loose their sense. 

The calculation procedures described above including ones for computation of momentum 
recoupling matrices and Talmi-Moshinsky-Smirnov coefficients were implemented in computer 
code. The ordinary procedures for generating of complete many particle sets for large spaces 
usually suffer of buildup of numerical errors. Their presence could even lead to undeffinite- 
ness of the obtained space dimension. Numerical instabilities do not concern the calculational 
procedures under preposition. All computations were performed with precise arithmetic using 
numbers represented in the form n/m/y/k, where n,m and k are integers. The application of 
the precise arithmetic is possible due to the orthogonal nature of the intrinsic density matrix 
formalism and the particulary character of developed calculational procedures that involve any 
numerical diagonalization and orthogonalization. The effectiveness of the proposed approach 
could be illustrated by the fact that intrinsic density matrices and related quantities of CFP 
nature for energy up to ^ and nuclei up to A = 11 could be calculated on 40 Mhz 

AT-386 within a few hours. 
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